Scaling Theory of the Mott Transition and Breakdown of the Griineisen Scaling Near 

a Finite- Temperature Critical End Point 



O 

(N 

C 
3 



Lorenz Bartosch, 1 Mariano de Souza, 2 and Michael Lang 2 

1 Institut fur Theoretische Physik, Goethe- Universitdt, 60438 Frankfurt am Main, Germany 
2 Physikalisches Institut, Goethe- Universitdt, 60438 Frankfurt am Main, Germany 

(Dated: April 27, 2010) 

We discuss a scaling theory of the lattice response in the vicinity of a finite-temperature critical 
end point. The thermal expansivity is shown to be more singular than the specific heat such that 
the Griineisen ratio diverges as the critical point is approached, except for its immediate vicinity. 
More generally, we express the thermal expansivity in terms of a scaling function which we explicitly 
evaluate for the two-dimensional fsing universality class. Recent thermal expansivity measurements 
on the layered organic conductor k>(BEDT-TTF)2 X close to the Mott transition are well described 
by our theory. 
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Scale invariance and universality are among the most 
striking concepts in physics [H-Q]. A prominent exam- 
ple is the critical end point of the liquid-gas transition of 
water at which the correlation length diverges and fluctu- 
ations occur on all length scales. Even though quantum 
mechanics can be important for understanding the na- 
ture of the different phases present, in the vicinity of a 
finite-temperature critical end point, a purely classical 
theory suffices to describe the transitions between these 
phases. As is well-known from the scaling theory of crit- 
ical phenomena, every critical end point is characterized 
by a universal scaling function and a set of critical expo- 
nents, known as the scaling exponents [iHl] . The scaling 
function as well as the scaling exponents do not depend 
on microscopic details of the underlying system. Instead, 
they are entirely determined by a few general properties 
like the dimensionality of the system, the symmetry prop- 
erties of the order parameter, or the presence or absence 
of long-ranged interactions. 

A very interesting class of materials consists of the so- 
called Mott insulators which have an odd number of con- 
duction electrons per unit cell and, according to band 
theory, should be metals. However, if the ratio of the 
on-site correlation energy U to the kinetic energy W be- 
comes larger than a critical value, all electrons localize. 
Prominent examples for the Mott transition are the Cr- 
doped V2O3 [H and the layered organic charge-transfer 
salts k-(BEDT-TTF) 2 A: jS^. Here, BEDT-TTF is 
bis(ethylenedithio)tetrathiafulvalene, hereafter abbrevi- 
ated ET, and X stands for various kinds of monovalent 
anions. As shown in Fig. [1] their ground states can 
be tuned by applying either chemical or external pres- 
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sure. While simple arguments based on symmetry 13 1 
suggest that the Mott transition lies in the Ising univer- 
sality class 14H17j|. this scenario is by no means obvious 
and is still under debate (l8l - [20| . Using the pressure- 
sweep technique to control the ratio U/W, conductiv- 
ity measurements by Limelette et al. on (Vi_ a: Cr x )203 
near the critical end point were reported to reveal mean- out that the conductivity can have a different scaling be- 
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FIG. 1. (Color online) Simplified phase diagram of the k- 
(ET) 2 X family as mapped out in Refs. _6-8]. While the salt 
with X = Cu[N(CN) 2 ]Cl (located at p = 0) is a Mott insu- 
lator, the hydrogenated variant of the X = Cu[N(CN)2]Br 
salt (ft8-Br) is a superconductor. Its deuterated variant (d8- 
Br), investigated in Ref. [l(J, is located at ambient pressure 
(dotted line) very close to the critical region. For simplicity, 
regions of inhomogeneous phase coexistence are left out. The 
thick, solid line corresponds to a first-order transition. It ter- 
minates in a critical end point (red dot), above which there 
is only a crossover between the Mott insulator and the metal 
(dashed line). The crossover line differs from the dashed- 
dotted line t w 1.74|ft| 8 ' 15 where signatures in the lattice 
response are expected to be the strongest. 



field behavior, with a crossover to scaling properties with 
3D Ising critical exponents However, unconventional 
critical behavior was claimed by Kagawa et al. for pres- 
surized quasi-2D k-(ET) 2 Cu[N(CN) 2 ]C1 (abbreviated k- 
Cl) [§| . The critical exponents obtained from the conduc- 
tivity measurements obeyed scaling relations but did not 
match any established universality class. In particular, 
Ising criticality was ruled out by the authors. 

A different explanation of the experiments by Kagawa 
et al. was given by Papanikolaou et al. 
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havior than the order parameter m itself. In particular, 
for the k-(ET)2X family, it was argued that the conduc- 
tivity E should scale like the energy density, implying 
<5E oc \m\ e , where 6 = (1 — ct)//3 is a combination of the 
critical exponents a and /3 related to the singular behav- 
ior of the specific heat and the order parameter, respec- 
tively With 9 = 8 for the 2D Ising universality class, it 
was shown that the conductivity measurements in Ref. Q 
are actually consistent with the assumption of Ising crit- 
icality. A complementary experiment, which concerns 
the anomalous response of the lattice in the vicinity of 
the Mott transition of deuterated K-(ET) 2 Cu[N(CN) 2 ]Br, 
was reported in Ref. [lfj. The data, which were ana- 
lyzed by assuming Griineisen scaling between the ther- 
mal expansivity a p (T) and the specific heat c p (T) [i.e., 
ct p {T) — T c p (T)\, were found to be at odds with the 
experiment by Kagawa et al. [9]. As the expansivity 
measurements were performed very close to but not ex- 
actly at the critical point observed in Ref. [i| , a definition 
of a critical exponent is problematic. In this Letter, we 
demonstrate the breakdown of Griineisen scaling in the 
vicinity of a finite-temperature critical end point. By 
elaborating on a scaling theory of the expansivity we re- 
analyze the expansivity data of k-(ET)2A and find con- 
sistency with 2D Ising criticality. 

By identifying the pressure p and the temperature T 
as independent parameters to fine-tune the phase transi- 
tion, we can already infer that the change of volume AV 
at the Mott transition is proportional to the order param- 
eter m. This is a simple consequence of the fact that the 
volume V and pressure p are conjugate variables, and it 
implies that the change in volume across the Mott tran- 
sition should scale as AV^ oc (T c — T)' 3 , where = 1/8 for 
the 2D Ising universality class. Apart from the different 
value of the critical exponent /?, this is quite analogous 
to the familiar liquid-gas transition. 

To set up a complete scaling theory, let us consider 
the singular part of the Gibbs free energy, f(t, h), which, 
for the 2D Ising universality class, can be written in the 
scaling form (see e.g., Ref. [211 ] ) 

f(t, h) = — In t 2 + \h\ d/Vh <5> (t/\h\ yt/yh ) . (1) 

Here, * - {T -T c - ({p- Pc ))/T Q and h = (p — p c — X(T — 
T c ))/p are temperature-like and pressure-like scaling 
variables which also contain small linear mixing terms. 
Apart from the critical temperature T c and pressure p c , 
the scaling variables depend on the four nonuniversal con- 
stants A, Ci and po- The finite values of A and £ 
account for the tilted scaling axes shown in Fig. [1] For 
example, A = dp/dT can be identified as the value of 
the finite slope at the end of the first-order transition 
line. Let us for now assume that in suitable units the 
parameters A and £ are sufficiently small such that it is 
reasonable to use dt/dp ~ and dh/dT » 0. As further 
explained below, this reasoning can be justified except 



for the nearest vicinity of the critical point. 

While the renormalization group eigenvalues yt = 1 
and yh = 15/8 for the 2D Ising universality class are 
known from Onsager's famous solution, the scaling func- 
tion Q(x) is not known analytically and has to be evalu- 
ated numerically. This, however, can be done with very 



high precision 21 



22] . All thermodynamic critical expo- 
nents can be derived by taking derivatives of Eq. (JlJ with 
respect to t or h. For example, if V, S, and G denote 
the volume, the entropy, and the Gibbs free energy of the 
system, respectively, the specific heat is given by 
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For A = h = this reduces to the well-known result for 
the singular part of the specific heat, 
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which, due to the absence of a power of t, implies a 
critical exponent a = 0. Of central importance here 
is the singular part of the thermal expansivity a p — 
V^ 1 (dV/dT) p = —V~ x (dS/dp) T , 
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which for h — > ±0 and t < scales as (see also Ref. [171 ]) 
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oc sgn(/i) (-t)- 1+ P . 



The generalized Griineisen ratio T P (T) is defined as 
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In a different situation, where there is a line of second- 
order transitions (but no critical end point), the pressure 
p enters only via the critical temperature T c (p), such 
that taking derivatives with respect to p is essentially 
the same as taking derivatives with respect to T. As a 
consequence, the Griineisen ratio stays constant even if 
c p (T) and cx p (T) diverge. However, this is quite different 
at a critical end point where t and h are independent pa- 
rameters. While it is, by now, well established that the 
Griineisen ratio can diverge at a quantum critical point 
23], a strong increase in T p (T) is also found for a critical 
end point at finite temperature, where for t < (and 
A = 0) 



"sing 



- sing 



oc sgn(/i) (-t) 



(7) 



h=±0 



It should also be noted that for h = and t > the 
singular part of the Griineisen ratio vanishes. A more 
complete picture of the vicinity of the critical point can 
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be derived from the scaling function. To this end, we 
combine Eqs. ([1]) and ((4]) to obtain the scaling form of 
the expansivity, 

a sing (t, h) oc sgn(h) \h\- 1+ ( d -y^y»y a (t/\h\y^). (8) 

Here, d is the dimension of the system. The scaling func- 
tion *!/ a (x) is given by 

$ (x) = _ B.a:*"(a:). (9) 

A plot of ^q(x) for the 2£) Ising universality class is 
shown in Fig. [2] [24|. \]/ a (x) assumes its maximum not 
at x = 0, but at the finite value x ~ 1.74, such that 
for fixed ft, the strongest lattice response is expected for 
t w 1 . 74 1 Tt, 1 8 / 15 ; which is indicated by the dashed-dotted 
line in Fig. [TJ Using the asymptotic behavior of ^ a (x) 
for large values of x|, 

*a(l) OC ( { Z% 7/8 > X _t (10) 
[ X ' , X — ?> CO, 

we recover Eq. ([5]) and a S in g = for h = and £ > 0, 
as expected. More generally, the singular part of the 
expansivity a S i ng can be calculated for arbitrary values 
of t and h, see Fig. G3 The strong lattice response at 
t w 1.74|/i| 8 / 15 and the sign change of the expansivity 
across the h = line can clearly be seen. 

Having established a general scaling theory for the sin- 
gular part of the thermal expansivity, we now analyze the 
thermal expansivity measurements of Ref . [loj . For sim- 
plicity, we assume that the nonsingular background con- 
tribution to the thermal expansivity can be approximated 
by a linear term, such that for the 2D Ising universality 
class we have 



a{T lP ) = Asgn[p - Pc - X(T - T c )} 
x [p- Pc -X(T-T c )]- 7/15 




FIG. 2. Plot of the scaling function ^(x) for the 2D Ising 
universality class. 




FIG. 3. (Color online) Plot of the thermal expansivity a s in g 
as a function of t and h which diverges at the critical end 
point t = h = 0. The line of first-order phase transitions 
along the negative t-axis [where the volume of the system 
changes according to AV oc (— t) 1 / 8 ] is clearly accompanied 
by a branch cut in the thermal expansivity. Note that due to 
the scale invariance near the critical point this plot looks the 
same on all length scales, i.e. rescaling t b vt t and h — ► b Vh h 
with arbitrary b we obtain exactly the same graph. 



For other universality classes, the exponents 7/15 and 
8/15 should be replaced by (l-/3)/(/3+7) and l/(/3 + 7 ), 
respectively, where 7 is the critical exponent associated 
with the compressibility k. A fit to the data published 
in Ref. [l(J by Eq. (pi} is shown in Fig. |4] The data 
are nicely described by our theory which demonstrates 
that the expansivity measurements of Ref. [lfj are in fact 
consistent with the assumption that the Mott transition 
in the k-(ET)2A family is in the 2D Ising universality 
class. 

Finally, let us explore the linear mixing terms in 
some more detail. For finite A, the specific heat c p 
and the thermal expansivity a p contain terms propor- 
tional to (XTq/pq) 2 k and (XTq/p )k, respectively, with 
k being more singular than c s j ng and a s j ng computed 
above. Consequently, even if (XTq/pq) is small, there 
is a region where this term becomes important and the 
Griineisen ratio saturates. A simple analysis shows that 
the crossover scales for the 2D Ising universality class are 
\T-T C \ » A 8 / 7 /S 15 / 7 and \p - p c \ w (X/B) 15 / 7 . How- 
ever, since for the experiment in Ref. 10] p — p c ~ 50 bar 
3> (A/-B) 15 / 7 s» 0.1 bar, having neglected the term pro- 
portional to Xk in our analysis above is clearly justified. 
It has also been argued that a finite shear modulus in 
solids can lead to a suppression of critical fluctuations 
near the critical point if the order parameter is propor- 
tional to AV (see Ref. [2(| and references therein). 

In summary, we have developed a scaling theory for 
describing the singular part of the thermodynamic ex- 
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FIG. 4. (Color online) Fit to the expansivity data for d8-Br 
crystal # 1 (diamonds 0) and # 3 (triangles V) from Ref. 
by Eq. (|lip . The fits turn out to be only weakly sensitive 
to the parameter A. A least-squares fit gives, for the fitting 
parameters of crystal # 1, T c + ((p—p c ) = 27.5 K, (p—p c )/A = 
26.7K, A/X 7/15 = 874 ■ lfr 6 K- 8/15 , B/\ 8/1B = 3.88K" 7/15 , 
C = 13.2 • 1CT 6 K _1 , and D = 0.85 • 1CT 6 K~ 2 . The straight 
dashed lines denote the linear background contribution. 



pansivity in the vicinity of a finite-temperature critical 
end point. As a salient result of our study, we find for 
h = ±0 and t < a universally diverging Griineisen ratio 
Tsing = "sing/csing oc sga(h) (—t) ~ 1+a +l 3 , except for the 
immediate vicinity \T-T C \ < (A/ B/ 3 +7)i/(/3+7-i) j w h e re 
r p (T) is expected to saturate. For the Mott transition in 
the vicinity of the critical end point, volume changes are 
a natural consequence of thermodynamics and are similar 
to the liquid-gas transition, implying AV oc (T c — T)' 3 . 
By explicitly evaluating the scaling function for the 2D 
Ising universality class, an excellent description of ther- 
mal expansion data for the layered organic conductor k- 
(BEDT-TTF) 2 A has been achieved, consistent with the 
conjecture made by the authors of Ref. [l7| and their 
analysis of the conductivity measurements in Ref. 
We stress that the approach presented here does not 
rely on any assumption regarding the system's under- 
lying microscopic interactions, including a possible in- 
volvement of lattice degrees of freedom [10( or system 
anisotropics. Our theory also predicts large lattice ef- 
fects at t ~ |/i| 8 / 15 , reaching out to both the low- and 
high-pressure sides of the critical end point. The obser- 
vations of a large anomaly in the thermal expansion of 
/i8-Br|25[, as compared to a tiny feature in the specific 
heat [261 ] , and a pronounced dip in the elastic constant of 
pressurized k-CI on both sides of the critical end point 
8] are consistent with the latter prediction. In fact, by 
extending our scaling theory to the singular part of the 
compressibility, filing oc d 2 f/dh 2 , we find that the dip in 
the sound velocity observed in Ref. Q should follow a 
line close to the dashed-dotted line in Fig. [TJ Note that 
these lines should not be confused with the crossover line 
obtained from the conductivity measurements of Ref. [t| . 
Finally, let us emphasize that due to its universal char- 



acter, our scaling theory should apply to a broad range 
of systems, including e.g. metamagnetic materials, for 
which an anomalously large Griineisen ratio and a sign 
change of the expansivity have recently been reported 
27] • 
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